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Abstract 

Symmetry-based ideas, such as the quark-lepton complementarity (QLC) principle and the tri- 
bimaximal mixing (TBM) scheme, have been proposed to explain the observed mixing pattern of 
neutrinos. We argue that such symmetry relations need to be imposed at a high scale A ~ 10 12 
GeV characterizing the large masses of right-handed neutrinos required to implement the seesaw 
mechanism. For nonhierarchical neutrinos, renormalisation group evolution down to a laboratory 
energy scale A ~ 10 3 GeV tends to radiatively break these symmetries at a significant level and 
spoil the mixing pattern predicted by them. However, for Majorana neutrinos, suitable constraints 
on the extra phases 02,3 enable the retention of those high scale mixing patterns at laboratory 
energies. We examine this issue within the Minimal Supersymmetric Standard Model (MSSM) 
and demonstrate the fact posited above for two versions of QLC and two versions of TBM. The 
appropriate constraints are worked out for all these four cases. Specifically, a preference for ot2 ~ vr 
(i.e. mi ~ —vri'i) emerges in each case. We also show how a future accurate measurement of #13 
may enable some discrimination among these four cases in spite of renormalization group evolution. 
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I. INTRODUCTION 

Outstanding recent experiments have increased our knowledge [l| of neutrino masses and 
mixing angles enormously. We are already certain that at least two of the three known 
neutrinos are massive, the heavier and the lighter of them being respectively > 0.05 eV 
and > 0.009 eV in mass. We also know that two of the three neutrino mixing angles are 
large: #23 ~ 45° and #12 ~ 34°, while the third is significantly smaller: # 13 < 12°. The 
total sum of the neutrino masses is also cosmologically bounded from above by 0(1) eV. 
Much remains to be known, though. The values of #13 and the leptonic CP violating Dirac 
phase Se, are still unknown. So is the ordering of the neutrino masses rrii (i = 1,2,3) - 
whether it is normal ({m^l > \mi^\) or inverted (Im^l < |m 12 |). We also do not know if 
the three neutrinos are hierarchically spaced in mass like charged fermions or if they are 
nonhierarchical. Our term nonhierarchical here includes both the inverted hierarchical (IH) 
case , i.e \m 3 \ <C |mi| ~ \m 2 \ ~ 0.05 eV, and the quasi-degenerate (QD) situation [3j, i.e. 
\mi\ ~ \m2\ ~ |m3| ^> 0.05 eV, the latter with either a normal or an inverted mass ordering. 
Neither of these scenarios is observationally excluded as yet and we focus on them. As per 
our present knowledge, the average neutrino mass could still in fact be anywhere between 
half of the atmospheric oscillation mass scale, i.e. ~ 0.025 eV and a third of the cosmological 
upper bound, i.e. ~ 0.3 eV. Finally, most theoretical ideas expect the three neutrinos to 
be Majorana particles whose masses m; can be complex. In that case, since one of their 
phases can be rotated away, there are two additional, possibly nonzero, phases [2] on which 
we do not have any direct information at present. This is because no convincing evidence 
exists as yet of neutrinoless nuclear double beta decay which is the only known direct probe 
{4] on these phases. Indirectly, of course, some constraints on these phases may also arise 
from considerations of leptogenesis Q]. It is nevertheless worthwhile to try to constrain 
these phases in some other way. That is one of the aims of the present work, which is 
an elaboration of our earlier shorter communication jy] with many additional results. In 
particular, we demonstrate here that, given the constraints on these Majorana phases, a 
measurement of #13 can make some discrimination among four scenarios considered by us 
despite renormalization group (RG) running. 

The observed bilarge pattern of neutrino mixing has led to the idea of some kind of a 
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symmetry at work. Several symmetry-based relations 1 have in fact been proposed, which give 



hat 



rise to specific neutrino mixing patterns. Two of the most promising mixing patterns, t 
we will be concerned with here, are (i) quark-lepton complementarity (QLC)JtI [g, [I], 1(3, 
and (ii) tribimaximal mixing (TBM) [12j. QLC involves bimaximal mixing [13] followed by 
the unitary transformation of quark mixing. A bimaximal mixing can in turn be generated 



symmetry 



by a /i-r exchange symmetry 14], an L M — L T gauge symmetry 15], or an S3 permutation 



The second step is inspired by SU(5) or SO(10) GUT, as discussed later. 



1* 



family 



Similarly, a tribimaximal mixing pattern may be obtained from an [17| or S3 
symmetry. However, a major issue in connection with such symmetries is the scale at which 
they are to be implemented. Neutrino masses and mixing angles are related directly to 
the corresponding Yukawa coupling strengths which run with the energy scale. There is 
as yet no universally accepted explanation of the origin of neutrino masses, but the seesaw 
m echan 1Sm Q , th e m0 , «eva bl e candidate so fa, The tan „ t tbe Ug nt nenMno 
mass matrix in family space in that case is M. v = —(m®) T 'M^mf?, where is the Dirac 
neutrino mass matrix (analogous to the charged fermion ones) and Mr the mass matrix 
for very heavy right chiral singlet neutrinos. If the Dirac mass of the heaviest neutrino is 
taken to be 1 - 100 GeV, the atmospheric neutrino data require typical eigenvalues of Mr 
to be in the 10 n -10 15 GeV range |20j. This is also the desirable magnitude for Mr from 



the standpoint of a successful leptogenesis 21] . From these considerations, we choose to 
implement the above mentioned symmetries at the scale A ~ 10 12 GeV. One can take issue 
with the particular value chosen for A. However, our conclusions are only logarithmically 
sensitive to the precise value of this scale. 

A question arises immediately on the application of such a high scale symmetry on the 
elements of the neutrino mass matrix. It concerns their radiative breaking via RG evolution 
down to a laboratory energy scale A ~ 10 3 GeV. The actual evolution 22|, |23j, |24J needs 



to be worked out in a specific theory which we choose to be the minimal supersymmetric 
standard model (MSSM 25j]). That is why we have taken A to be of the order of the explicit 
supersymmetry breaking or the intra-supermultiplet splitting scale O(TeV). Once again, our 
calculations are only logarithmically sensitive to this exact choice. The point, however, is 



1 Here one should perhaps make a distinction between a symmetry of the Lagrangian and just a special rela- 
tion among coupling strengths or masses. Nevertheless, the relations of concern to us can be implemented 
through specific symmetries of the Lagrangian. 
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that - for nonhierarchical neutrinos - symmetry relations formulated at A will in general 
get spoilt on evolution down to A. 

The full RG equations for the evolution of neutrino masses and mixing angles in the MSSM 
have been worked out 0, 



24| in detail. In particular, the evolution effects on the mixing 
angles are found to be controlled by the quantities [6| A r |mj + mj| 2 /(|mj| 2 — |mj| 2 ) where A T , 
to be specified later, is a small fraction <10~ 2 , while i,j refer to the concerned neutrino mass 
eigenstates. Consequently, these effects are negligible for a normal hierarchical mass pattern 
with 1 777.3 1 ^ \ m 2\ ^ \ m i\- RG effects can become significantly large only when neutrinos are 
nonhierarchical. There is another important characteristic of the above-mentioned ratios. 
While their denominators involve only the absolute masses the numerators involve the 
combinations \rrti + rrij\ 2 . Therefore, with appropriate constraints on the neutrino Majorana 
phases, the desired symmetry relations can be approximately preserved at the laboratory 
scale A even for nonhierarchical neutrinos - in agreement with the mixing pattern that 
has emerged from the oscillation data. The constraints on the majorana phases and the 
consequent discrimination among the scenaios by a measurement of #13 constitute our main 
results. Our work is somewhat complementary to that of Ref. 2(| in the QLC sector and 
Ref. [271 in the TBM sector. 

In this paper we work out in detail the last-mentioned constraints on the neutrino Majo- 
rana phases in the (i) bimaximal mixing + QLC and (ii) TBM scenarios respectively. Each 
of these comes in two variations. So we have in all four cases at hand. Thus, the scope of 

n 

the present work is much larger than our earlier shorter communication [6] which addressed 
only one version of QLC and did not consider the implications for #13. A major technical 
observation utilized by us is the following. Suppose #f 3 , the high scale value of the angle 6> 13 , 
is sufficiently small (as is the case for the situations considered here) such that (9 (#13) terms 
can be neglected in comparison with other 0(1) terms in the RG equations [24] . Then the 
neutrino mass matrix at the laboratory scale A, becomes analytically tractable in terms 
of its high scale form A4^. In fact, the relation obtained looks quite simple and transparent. 
The step from there to explicit constraints on the neutrino Majorana phases is then shown 
to be quite straightforward. The rest of this paper is organised as follows. Sec. [Til contains a 
description of the parametrisation that we find convenient to adopt for nonhierarchical neu- 
trino masses. In Sec. IIII[ we introduce two versions each of the QLC and TBM scenarios to 
be implemented at the high scale. In sec. IIVI we discuss the energywise downward evolution 
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of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix in general, and its effects 
on the predictions of the scenarios under consideration. In sec. |V]we study the constraints 
on neutrino and MSSM parameters in order for the scanarios to be valid and explore if these 
scenarios may be distinguished by means of more accurate measurements of the neutrino 
mixing angles. The concluding sec. IVII consists of a summary and the discussion of our main 
results. 



II. PARAMETRISATION OF NONHIERARCHICAL NEUTRINO MASSES 



We work in the convention [28] in which the neutrino mass eigenstates |z/ 2 ), \v 3 ) are 
related to the flavour eigenstates \u e ), \u T ) with the unitary mixing matrix U u : 



I fa) = U ai \Vi) 



(1) 



a and i being flavour and mass indices respectively. We take the neutrino mass term in the 
Lagrangian to be 



1- 



Thus, 



' m 1 » 



(2) 



utM v u: 



v 



(3) 



m 2 
m 3 

where are in general complex. However, one of the three phases of 7711,2,3 can be absorbed 
in the overall phase choice of in (T5]). We can therefore choose 



mi — \ m i\ j ^2 — |?^2|e ia2 , m 3 = \m 3 \e za3 



where «2,3 are real. Experiments with atmospheric neutrinos tell us that 2 81 129] 



\m 3 \ 2 - \m 2A \ 2 \ = (2.4 ± 0.3) x 10~ 3 eV 2 , 



(4) 



(5) 



while experiments with solar electron neutrinos and reactor electron antineutrinos yield 
28, £] 

5m 2 s = \m 2 \ 2 - |mi| 2 = (7.9 ± 0.4) x 10~ 5 eV 2 . (6) 
For charged fermions (/ = u, d, I), the mass term is 



1- 



-^rnass - -^fRamfapfip + fl.C. 



(7) 



The corresponding mass matrix m/ is put into a diagonal form by 



Ujm\rrifUf = \m K f 



P)|2 



(8) 



Now the unitary Cabbibo-Kobayashi-Maskawa (CKM) and the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrices, whose elements contribute to the observed quark and 
neutrino processes respectively, are given by 



CKM 



ulu. 



f^PMNS — UlU, 



(9) 



One can write £7pmns m the standard basis [28( in terms of the angles #12, #23, #13 and the 
CP violating phase 5e- With = sin% and c^- = cos%, 



/ 



U 



PMNS 



C12C13 



S12C13 



-S12C23 - Ci2S 2 3Sl 3 e^ C12C23 - s 12 s 2 3Si 3 e l5)! 



\ 



S12S23 - Ci 2 C23Sl 3 e 



iSr 



-C12S23 - Sl2C23Sl 3 e 



si 3 e- i5e 

S23C13 
C2 3 Ci 3 



(10) 



The experiments mentioned earlier then also tell us that 28|, [29] 



012 = 33.9°±1.6 C 
e 23 = 43.3°±8.2 C 
0n < 12° . 



(11) 



We find it convenient to parametrise the absolute masses \rrii\ for nonhierarchical neutri- 
nos in terms of three real parameters m ,pA and e$ as follows: 

I mi I = m (l - Pa)(1 - e s ) , 
\m 2 \ = m (l - Pa)(1 + es) , 

|m 3 | = m (l + pa) • (12) 

In eqs. f lT2|) . m defines the overall mass scale of the neutrinos, whereas pA and €5 are 
dimensionless fractions with — 1 < pa <C 1 and < eg < \pa\ for nonhierarchical neutrinos. 
The sign of pa is positive (negative) for a normal (inverted) ordering of neutrino masses. 
Moreover, pa ~ — 1 (|/0a| ^ 1) f° r the IH (QD) case; in either case e$ <C 1. For comparison, 
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FIG. 1: The parameters p A , es and T as well as the value of Y2 \irii\ (eV) as functions of mo, for a 
normal (left panel) and an inverted (right panel) mass ordering of neutrinos. Note that pa and V 
are negative for inverted mass ordering. For normal ordering, es goes to near unity with low mo- 
it may be noted that for normally hierarchical neutrinos, pa ~ es ~ 1- We can further write 
the solar and atmospheric neutrino mass squared differences as 

| 2 - |mi| 2 ^4m 2 (l ~ p A ) 2 t s , 

|2 /I™ I/O 1 I™ I 



Sm 2 s 



\5m 2 A \ 



4m \p A \ 



|m 3 | 2 -(|m 1 |/2+|m 2 |/2)^ 
Utilizing ([6]), ([5]) and (|T3|) . we see that 

m > 0.024 eV . 

Also, the cosmologically bounded sum of neutrino absolute masses is given by 

£j|mj| = 3 mo(l — Pa/3) < 1 eV . 

From ([13]) and ([T5D, it follows that 

4(EKI) 2 _ (1-Pa/3) 2 
\Pa\ 



(13) 



(14) 



(15) 



(16) 



9 \5m 2 A \ 

Utilising §5§, (TTB"j) and the cosmological upper bound (TT5"|) . we get \p A \ > 5.5 x 10~ 3 . 

In FiglU we show how pa, es and Yl l m «l behave as functions of m . We also find it 
convenient to define the derived dimensionless parameter 



T = — 

Pa 



Pa 



(17) 



whose behaviour is included in the figure. Since \pa\ < 1, the sign of T is positive (pa > 0) 
for a normal ordering and negative (pA < 0) for an inverted ordering of the neutrino masses. 
The inequality 5.5 x < \pa\ < 1 translates to |T| < 180. We can generate typical numbers 
for nonhierarchical neutrinos: In the IH scenario, for instance, if mo is chosen as 0.025 eV, 
we get pa ^ — 1, T « 0~, and es ~ 8 x 10~ 3 . On the other hand, for QD neutrinos, choosing 
m = 0.2 eV, we have \pa\ ~ 1-5 x 10~ 2 , \T\ ps 65, and eg « 5 x 10" 4 . 



III. SYMMETRIES AND MIXING ANGLES AT THE HIGH SCALE 

The neutrino mass matrix at the high scale A originates from a dimension-5 operator 

In f fT8l) . £ a and H are the SU(2) doublet lepton and Higgs fields respectively and c a p are 
dimensionless coefficients that run with the energy scale. Then 

Mt aP ~ c ap V - , (19) 

where v = 246 GeV and A ~ Mmaj, the Majorana mass characterising the heavy SM-singlet 
Majorana neutrino N. The symmetries at the high scale A give rise to specific structures for 
the matrix A4^, and hence predict the values of the mixing angles 9^ at the scale A. There 
is an issue about a consistent definition of v\p at all scales. Solar neutrino experiments tell 
us that jm^l > where A is the laboratory scale. We define v\ and v<i at higher scales 

in a way such that |m 2 | > |mx| at all scales, and in particular, \rri2 \ > |m^|. 

In this section, we introduce four different symmetries and the corresponding predictions 
on the neutrino mixing angles at this high scale: 



A. QLC1 



Quark-lepton complementarity 



ll) links the difference between the measured 



and the maximal (Le. 45°) values of the neutrino mixing angle 812 to the Cabbibo angle 



r 



12.6° ± 0.1 c 



281 ] . We first follow a particular basis independent formulation "QLC1" 
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9J of this principle 2 : 

f^PMNS = ^CKM^bm , (20) 

where U v ^ m is the specific bimaximal form [l3] for the unitary neutrino mixing matrix. 
Eq. (JH]) gives rise to the "QLC1" relation 3 

0t 2 + ^ = \ + O{0l) . (21) 

The identification of (1201) as a statement of QLC becomes more transparent in the basis 
with U u = I, i.e. where the matrix Y^Y U is diagonal. It follows from (jHJ) that Vckm = Ud 
in this basis. Now a comparison of (Q and (1201) . together with the assumption of U v being 
Uvfrxai yields the SU(5) GUT-inspired quark-lepton symmetry relation Ud = U\. Eq. (|20|) . as 
it stands, is basis independent, however. 

Eq. (|20l) yields the neutrino mixing angles at the high scale A to be 

4 = \ - A + 0(61) « 35.4° 

# 2 A 3 = |-|K*|-f + C(^ 3 )~42.1°, 

0f 3 = A + 0(fl c 3 ) « 8.9° . (22) 
Thus, QLC1 predicts a value of 0f 3 that is close to the current experimental bound. 



B. QLC2 

In a second version of quark-lepton complementarity, "QLC2" {§], [^] , one assumes a 
bimaximal structure for the charged lepton mixing matrix, Ue = U^ m and the form 

t^PMNS = Ue t b m Vc KM , (23) 

for the PMNS matrix. Eq. ( 1231) yields in a straightforward way the relation 

0&+0c=1 + 0{9l) . (24) 

2 There could be a more general statement of QLC1 with an additional diagonal phase matrix T$ between 
Vq KM and U Vt bm- But for consistency and simplicity, we choose eq. (j20|) . 

3 We do not distinguish between 9^ and 9* since the running of 9 C is negligible on account of the hierarchical 
nature of quarks belonging to different generations. 
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One may note that, in the basis with Ud = I, i.e. where YjYj is diagonal, fl23|) yields the 
SO (10) GUT-inspired relation U u = U u . 

Eq. f[2"3"j) leads to the following values for the nautrino mixing angles at the high scale: 



#12 - 


7T 

4 ~ 


9 C + O(0 3 c ) « 32.4° , 




#23 = 


71 

4 ~ 


+ C(^) « 43.4° , 

v2 




y i3 - 


\Vcb\ 

s/2 


+ O{0l) « 1.6° . 


(25) 



The value of 0f 3 predicted in QLC2 is beyond the measuring capacity of the neutrino exper- 
iments planned during the next decade. 



C. TBM1 



U A = — 

i/,tbm /g 



^(^^^^(sin- 1 -^). 



(26) 



The tribimaximal form of the neutrino mixing matrix is given by 

^2^0^ 

i 

-1 V2 
K 1 -V2 V3 j 

In the standard TBM scenario [l2j, which we refer to as TBM1, one has Up MNS = U„ 
since the charged lepton mass matrix at the high scale is already flavour diagonal. Then we 
have 

0° . (27) 



a 12 ~ 



35.3° , 0* = 45° , 9[\ 



D. TBM2 



Small deviations from the tribimaximal scenario TBM1 above have been considered in 
the literature 



34 



35j , where the deviation originates from the mixing in the charged lepton 



sector. Here we consider the version in Ref. 
where Vn has the form 11] 



/ 



V, 



IL 



i e c /3 o 

9J3 1 -\V cb 
\V cb \ 1 



34], and call it TBM2. Here C/pmns — VaJ^^tbm) 
\ 



+ 0(6, 



c) ' 



(28) 



/ 



10 



with the factor of 1/3 coming from the Georgi-Jarlskog relation 36( m M /m s = 3 at the GUT 
scale. As a result, we have at the high scale 



9 



12 



#23 



#13 



sin" 1 -L - + 0(6l) « 32.3° 
>/3 3^ 



7T 



--|U c6 | + 0(# 3 ) ^42.7° , 



0, 



3^ 



0(£ 3 



3.1° 



(29) 



IV. HIGH SCALE STRUCTURE AND DOWNWARD EVOLUTION 



As explained in the Introduction, our idea is to start with a specific structure of the 
neutrino mass matrix Ai^ that is dictated by some symmetry at a high scale A ~ 10 12 GeV. 
We would then like to evolve the elements of Ai u down to a laboratory energy scale A ~ 10 3 
GeV. This involves studying the (one-loop) RG evolution of the coefficient functions c a p in 
(j!8p between A ~ 10 12 GeV and A ~ 10 3 GeV. In case the considered high scale neutrino 
symmetries are consequences of grand unification, we need to assume that the threshold 
effects 30[ between the GUT scale ~ 2 x 10 16 GeV and A are flavor blind so that they do 
not spoil the assumed symmetry relations in the downward evolution from Mgut to A. We 
also note that effects of evolution on the masses and mixing angles of charged fermions are 



known 



3jj to be negligibly small 4 on account of the hierarchical nature of their mass values. 



At one loop, the neutrino mass matrices at the scales A and A are homogeneously related 



23, 



33|]: 



m: 



I K Xl Mt 1 K 



where 



exp 



t(A) 



K{t)dt 



t(A) 



(30) 



(31) 



is a scalar factor common to all elements of M.„. In eq. (1311) . t(Q) = (167r 2 ) 1 hx(Q/Qo) 
with Q (Q ) being a running (fixed) scale, and the integrand is given by 



K{t) = -6g 2 2 (t) - 2g 2 Y {t) + 6Tr (Y^Y u ){t) 



(32) 



4 The value of \V c b\ does run by about 0.01 in the MSSM due to the top quark U(l) coupling. However, at 
the level of accuracy that we are concerned with, this is inconsequential. 
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in a transparent notation, g 2 y being the SU(2)l, U(l)y gauge coupling strength and Y u the 
up-type Yukawa coupling matrix. Finally, the matrix I K has the form 



I K = exp 



t(A) 



t(A) 



(33) 



Yg being the Yukawa coupling matrix for charged leptons. 

Although some of the neutrino mixing matrices in various scenarios in Sec. IHII have been 
motivated in terms of grand unification in bases where the symmetries involved may be 
clearly observed, for the RG evolution of all scenarios we choose to work in the basis where 
the charged lepton mass matrix is diagonal. In this basis, 

Y}Y t = Diag (yl y% y 2 T ) . 

We can neglect y 2 in comparison with y 2 in (l3~4l) to get the result 

Z K « Diag(l, 1, e-^) = Diag(l, 1, 1 - A T ) + 0(A^ 



(34) 



4), 



where 



32 ] 



t(A) 



\y T (t)\ 2 = m^(tan 2 /3 + l)(87rV)- 1 ln(A/A) . 



(35) 



(36) 



t(A) 



Here v = \jv\-\-v\ and tan/3 = v u /vd, where v u (vd) is \/2 times the vev of the up (down) 
type neutral Higgs scalars. Evidently, A r is a small number for the allowed range of tan/?: 
e.g. A T pa 6 x 1CT 3 for tan/3 = 30, justifying our neglect of the 0(A 2 ) terms 5 . 
The substitution of fl35l) into fl30l) leads us to 



M v oc 



/ 



V 



1 
1 
1 — A T 



\ ( 



J 



1 
1 
1 - A T 



/ 



+ 0(A 2 T ) 



(37) 



where the proportionality is through the scalar factor Ik given in ( !3~TI) . It must be emphasised 
that (1371) is valid only when Ai v is written in the basis where the charged lepton mass 
matrix is diagonal. 

The matrix Ai„ is complex symmetric. Writing it in the general form 



Mi 



( A B 
B D E 
C E F 



(38) 



Note that a mistake of a factor of 2 in eq. 9 of [6| has been corrected here. 
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we have 



( A B C(l - A T ) ^ 



+ 0(A 2 T ) . (39) 



Mt = B D E(l — A T 

yC{l-A T ) E(1-A T ) F(1-2A T ) J 

Since both M A,X are complex symmetric matrices, they can be diagonalized as 



(U A y M A (U A )* = Diag(K|,K|,K|) 
(U x y M X (U X )* = Diag(K|, Kl, Kl) • (40) 

Note that, since we are working in a basis where the charged lepton mass matrix is diagonal, 
jjA (A) j g same as the net leptonic mixing matrix C/pmns ® at the scale A (A). 
The unitary matrix £/pmns m ay be parametrised in its most general form as 

[/pmns = Diag(e^,e^,e^)^23(^3)Diag(l,l,e^)^i3(^i3) x 

Diag(l, 1, e-^)^i 2 (^i 2 )Diag(l, e~^\ e^ 2 ) , (41) 

where Rij(6ij) is the matrix for rotation through the angle dij in the i — j plane, 5? is the 
CP violating Dirac phase, 02,3 the Majorana phases and 4> e ,fi,r the so-called "unphysical" 
additional phases required to diagonalise the neutrino mass matrix. Here we have already 
used the freedom of choosing a\ — 0. 

If 9 A 3 vanishes on account of the symmetry requirement at the scale A, the evolution of 
the mixing angles can be computed analytically in a simple manner. The matrix U A can 
then be written as 

U A = i? 2 3(^ A 3 )i?i 2 (^ 2 )Diag(l, e"^/2 ; ? (42) 

since the Dirac phase contribution vanishes and the phases </> e ,/i,r can anyway be absorbed 
in the charged lepton phases. RG evolution will modify the angles #23, #12 as well as the 
phases 0:2,3, at the same time generating nonzero values for the mixing angle #13 and the 
Dirac phase 5#. The phases <f> e iitT that may get generated can always be absorbed in the 
phases of the charged lepton flavour eigenstates. 

If we now approximate the deviation of U x from U A by the retention of only terms that 
are linear in A r , we can write the modified mixing angles as 

e x 2 = e A 2 + k l2 A T + o(a 2 t ) , e x 3 = e A , + k 23 A T + o(a 2 t ) , e x 3 = k 13 A T + o(A 2 T ) . (43) 
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The modified phases may similarly be written as 

«2,s = «2,s + «2, 3 A T + 0(A 2 T ) , 5 £ A = d e A T + 0(A 2 T ) , (44) 

where we expect 02,3 and di to be 0(1) quantities. In this paper we shall only be concerned 
about the deviation kijA T of the mixing angles Oij from their high scale values. We shall see 
a posteriori that |fcjjA T | is always much less than unity so that we can ignore its quadratic 
and higher powers. 

The values of kij are found to be 

k 12 = ^sin2^ 2 sin 2 ^ 3 





mf + m 2 


2 


\ m 2 


2 _ 


mf 


2 



k 23 = 1 sin 2^3 cos 2 e ^ZlV 1L + siu 2 ef 2 
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h 3 = - sm 26f 2 sin 20* 1 2 3 ' - 1 A 3 ' , 45 

4 \l m 3T — l m 2 I l m 3T — l m ll / 

where mf are the masses \mf \e ia ^ at the high scale. The RG evolution of \mi\ may be 
parametrised by 

I K \mf\ (l + niA T + C(A 2 )) (46) 



where I# is the scalar factor given in (131]) and are 0(1) numbers 24|. Then, taking m ; 
to be the masses \m^\e ia i at the low scale introduces an error of 0(A T ) in ky, and hence 
of O (A 2 ) in 0y. Eqs. ( 1451) are therefore valid even with the low scale values of the masses 
and Majorana phases. The same argument is true even for the mixing angles 0y. Since 
all observations in the neutrino experiments are made at laboratory energies, we henceforth 
drop the superscript A for all the quantities at the low scale. 

The expressions (l4"5j) h ave been derived starting with 9f 3 = 0, and they agree with the 
general expressions in [24j in the limit 0f 3 — > 0. For finite 0f 3 , the k^ given in (143!) have an 
error of O(0f 3 ), i.e. the mixing angles at the low scale have an error of O(0f 3 A T ), which 
we neglect in our analytical approximations. Since for all the scenarios considered here, 
9f 3 A T < 0.1°, this is a justified assumption. Thus, even with small nonzero values of 9f 3 , 
as is the case with QLC1, QLC2 and TBM2 scenarios, we can justifiably use the result 
6ij ps Of- + kijA T with k^ given by (145]) . Another requirement for the validity of eqs. (l4"5j) is 
that the values of and \nii\ 2 — \nij\ 2 should be described accurately by the 0(A T ) terms 
in their RG evolution. This condition does not impose any restriction on mf. However, 
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for A T > (|m 2 | 2 — |m A | 2 )/(m A ) 2 , the 0(A 2 ) terms dominate over the 0(A T ) terms in the 
1 777.2 1 2 — |^i 1 2 evolution [24|, which results in eqs. (|4"5l) breaking down. Thus for the validity 
of these equations, we require 

K) 2 -A T <K| 2 -K| 2 . (47) 

Eq. (|4T|) may be violated if \m 2 | 2 — |m A | 2 is indeed very small. However, even in such 
a situation, many of the qualitative features following from our analytical treatment will 
continue to remain valid. This will be confirmed by the exact numerical analysis given later. 

One can make some further observations on eqs. (I4"5l) . Since k 12 is always positive, 
the measured value of 6*12 can never be smaller than # A 2 . Moreover, in order to have any 
significant effect of RG evolution, some are required to be >O(10) since A r < 0(1O~ 2 ). 
This requirement is not met in the case of a normal hierarchy which has all \ky\ ~ and 
hence there is no significant RG effect there. Specifically, for an enhancement of \ki 2 \, |m A | 2 
must necessarily be ^> \m 2 | 2 — |m A | 2 . This inequality is obeyed whenever neutrinos are 
nonhierarchical, i.e. either in the IH or the QD case. On the other hand, the enhancement 
of both |fc 23 | and \ki 3 \ requires the inequality \p\\ <C 1, satisfied only by QD neutrinos. 
However, these are necessary conditions. The occurence (or not) of an actual enhancement 
depends on the Majorana phases a 2i3 too. In order to see this dependence explicitly, let us 
rewrite (T4"5T) in terms of these phases and our parameters p\ and of ffT2l . defined at the 
high scale. We then have, with finite 6>f 3 , 

k 12 = -^sin2^ 2 sin 2 ^ A 3 [l + cosa 2 v + (e^) 2 (l-cos^)]+0(^3), 

r A 

k 23 = — sin 26> A 3 [l + cos 2 6> A cos(a 2 - a A ) + sin 2 #f 2 cos a A ] 

+ ^sin2^ 2 sin20 2 A 3 + O(ei^3), 
k 13 = ^ sin 2# A sin 2fl A [cos(a A - a A ) - cos a A ] + 0(e A , A ) . (48) 

o 

It is clear from (T4H1) that the values of the Majorana phases a 23 are crucial in controlling 
whether the evolution terms kijA T are dangerously large or not. For instance, as a 2 — > 0, 
we have ki 2 « (4e^)~ 1 . For nonhierarchical neutrinos, the latter is likely to destroy any 
high scale symmetry statement on Qxi ■ When a 2 increases from zero, |fci2| decreases rapidly, 
becoming as small as le^l/4 when a 2 = n, i.e. m A « — m 2 . Moreover, fj48l also indicates 
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that both I&23I and | A;i3 1 are enhanced when either a 3 = 0, or a 2 = for a nonzero a 2 , 
though |fc 13 | gets highly suppressed when = 0. 

As noted in Sec. [TT1 the sign of T A is positive for a normal ordering and negative for an 
inverted ordering of the neutrino masses. Thus k 23 is always positive (negative) for a normal 
(inverted) ordering. However, the sign of A43 is controlled not only by sgn(r A ) but also by 
a combination of Majorana phases. 

In the remaining part of this section, we enumerate the predictions for the four sce- 
narios considered in this paper. On evolution down to the scale A, the mixing an- 
gles of the corresponding mixing matrix C/pmns are given in all the four cases by 8ij = 
Of, + kijA T + 0(9f 3 A T , A 2 T ) where the values of % are given by pgj) . 



A. QLC1 

On evolving ( 1221) to the laboratory scale, the net leptonic mixing angles are found to be 
0u = \ ~ ^ + k^ cl A T + <D{61 9t 3 A T , Al) , 



02 3 = '- - \V cb \ - I + 4 LC1 A r + O(0l ef 3 A T , A Tj , 







13 



°c 7.QLC1 A 

71 13 T 



+ o(e s c ,et 3 A T ,A 2 T ), (49) 



where the neglected terms are < 0.1°. The absolute value taken for the RHS of #13 is in 
order to keep to the convention of defining % > 



Since tt/4 - 6> c /x/2 w 35.4° and £^ 2 LC1 A T > 0, it follows from (ggj) that 6 12 in this scenario 
cannot be less than 35.4°. Further, since vr/4 - \V cb \ - 6 2 J4 w 42.1°, and A;^ LC1 A r > (< 0) 
for a normal (inverted) ordering of neutrino masses, a consequence of ( 1491) is that 623 is 
greater than (less than) 42.1° for a normal (inverted) ordering. Finally, the predicted value 
of 0i3 is 9 C /V2 w 8.9° in the absence of RG running, but it can be greater or less than 8.9° 
depending on the values of the Majorana phases. The detailed numerical analysis will be 
presented in Sec. [V] 
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B. QLC2 



After RG evolution, the high scale angles (1251) in this scenario evolve to 

012 = \ ~ Oc + k^° 2 A T - +0(6l et 3 A-r, At) , 

023 = t - ^ + 4 LC2 At + O(0 C 3 , 0f 3 A T , A?) , 



4 ^2 



9 



13 



it; 



cfe| 7.QLC2 



fi, 13 AA T 



+ 0{9* e ,e&A T ,&i), (50) 



v/2 

where the neglected terms are < 0.1° as before. 

The lower bound on # 12 in this scanario being 7r/4—9 c m 32.4°, significantly lower values of 
9 12 than QLC1 are allowed. Also, here 9 2 3 is greater than (less than) 7r/4— |V^,|/v2 ps 43.4° 
for normal (inverted) neutrino mass ordering. The major difference from QLC1 is in #13: 
the value #13 in QLC2 is only \V c b\/^/2 ps 1.6° in the absence of RG running. It can increase 
with RG running, but the extent of this increase is restricted from the observed 6* 12 values 
which restrict the values of the Majorana phases in turn. The detailed numerical analysis 
will again be presented in Sec. [V] 



C. TBM1 



The mixing angles at the low scale here are simply given by 

012 = sin- 1 -L + k™A T + ^ ^ ^ 
023 = \ + k™^A T + 0{9l9^Al), 

9 13 = A™A T + (9l 9t 3 A T , A r ) . (51) 

On similar lines to the arguments given for the QLC scenarios, here (i) the minimum value 
of 6*12 is sin~ 1 (l/v / 3) ~ 35.3°, (ii) the value of #23 is greater than (less than) 45° for normal 
(inverted) hierarchy, and the value of #13 vanishes in the absence of RG running. Since 
(|48p shows that # 13 does not run if the Majorana phases vanish, any observed deviation 
of 0i3 from zero in this scheme will indicate nonvanishing Majorana phases. The detailed 
numerical analysis appears in Sec. [V] 
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D. TBM2 



The mixing angles (1291) for this scenario evolve to the laboratory scale and become 
*i2 = sin- 1 i= - A= + k™A T + O{0l 9t 3 A T , A 2 ) , 



023 = \ ~ \V*\ + k™A T + O(0l 9t 3 A T , A^ 



tJ ' 



9 



13 



@c 7.TBM2 a 

3^ 13 T 



+ 0(9 3 c ,9t 3 A T ,A 2 T ). (52) 



The minimum allowed value of 9 i2 in this scenario is sin _1 (l/\/3) — 9 c /(3y2) ~ 32.3°, since 
/jTBM2 > q ^he V£ j ue f # 23 j s greater than (less than) tt/A — \V c b\ ~ 42.7° for normal 
(inverted) hierarchy. Finally, the value of # 13 is c /(3\/2) ~ 3.1° in the absence of RG 
running. Again, the detailed numerical analysis with RG running is presented in Sec. |V] 



V. CONSTRAINTS ON m , tan/3 AND MAJORANA PHASES 

In this section, we explore the current limits on the parameters of the four scenarios 
considered above. The main parameters governing RG running are mo, tan/3, the Majorana 
phases «2,3, and (to a smaller extent) the Dirac phase 8g. Our aim is to find the range 
of values of these parameters allowed by the current data. The four scenarios also lead to 
slightly different predictions for the mixing angles; and accurate measurements of these an- 
gles should distinguish among them in the absence of RG running. However, since the latter 
spoils high scale symmetries in general, and since the values of all the relevant parameters 
are still not known, the low scale predictions of the mixing angles are expected to ovelap. 
We explore in detail whether this still allows one to discriminate among the scenarios con- 
sidered in this paper. Moreover, we study the correlations between the deviations of the 
mixing angles from their high scale values. 

We use the 3a ranges for the neutrino mass and mixing parameters 

7 x 1(T 5 eV 2 < 5m% < 9.1 x 1(T 5 eV 2 , 1.7 x 10~ 3 eV 2 < \5m\\ < 3.3 x 1(T 3 eV 2 , (53) 

30° < 9 12 < 39.2° , 35.5° < 9 23 < 55.5° , 9 13 < 12° (54) 

at the low scale. At the high scale, we start with the values of 9^ dictated by the scenario 
under consideration, and a range of 5ml, A values that are consistent, after RG evolution, 
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with the low scale measurements (153]) and (1541) . We show our results for a normal mass 
ordering of neutrinos. However, the constraints in the case of an inverted mass ordering are 
almost identical. 

We give our constraints in terms of the values of m A and a A at the high scale. These 
are related to the low scale values of m and a 2 as indicated in fH6T) and (jUJ) respectively. 
In particular, the evolution of mo is controlled mainly by Ik ~ 0.71, which makes itlq m 
0.71 m A . The bounds on m A shown in the figures in this section can then be easily translated 
into bounds on the low scale value of mo- The RG evolution of a 2 is given by a\ ~ a 2 v + a 2 A T 
with [2J] 

^ « I 'a^T A.2 cos < sin 2 ^ 3 sin a A , (55) 
|m 2 1 — \m-y \ 

which may be used to translate the a A constraints to low scale values of oio- 

n 

We have neglected possible Planck scale effects [37J which may change the value of 6\2 by a 
few degrees for quasidegenerate neutrinos, leaving the other two angles virtually unaffected. 
Inclusion of these effects would relax [38] the constraints in the Wq - plane by a small 
amount. 



A. Limits on m A , a 2 and tan/3 

Out of the three leptonic mixing angles, #i 2 is the one measured with the greatest ac- 
curacy currently. Since all the scenarios have specific predictions for 612 in the absence of 
RG evolution, the measured value of 6* 12 can put the strongest constraints on the running 
parameters. From (H31 and ( 1481) . the running of Q\2 is expected to be independent of a A . 

We show in Fig. [2] the 3<r allowed regions in the m A — plane for two tan (3 values. 
This figure has been obtained from numerical solutions of the RG equations. The figure 
agrees very well with our expectations from the analytic expressions ( HHj) . On account of the 
occurence of the small quantity in the denominator of fci 2 in eq. fHHl) . strong constraints 
ensue on m A and a 2 . At large tan/3 (equivalent to a relatively large A r ), the value of a 2 
has to be near 7r (i.e. m A w — m A ) underscoring the necessity for a nontrivial Majorana 
phase. However, the requirement is less severe for a smaller tan f3. The coefficient a 2 in f l55l 



characterising the evolution of a 2 vanishes when = 7r [24j . As a result, the preferred 



value of a 2 , viz. = 7T, is equivalent to a 2 = 7r at all scales. 

The constraints displayed in the figure are calculated for two fixed tan/3 values. The 
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FIG. 2: Constraints from the 3<r allowed range of 612 on the m,Q (eV) - ot2 (radians) parameter 
space for tan (5 = 5 and 20 in the QLC1 scenario. The regions above the contours are excluded by 
data for that particular value of tan (3. The peak at ai ~ tt is noteworthy. 

results for other values of tan f3 may be extrapolated from the figure. However, we observe 
numerically that the 9 12 constraints on the - plane (Fig. [2]) depend essentially on 
the combination tan /3. Therefore, in Fig. [3j we show the allowed region in the tuq tan (5 
- a.2 parameter space for all four scenarios. The differences among the four scenarios arise 
primarily from the differences in the values of 9^ 2 . These are 35.4° (QLC1), 32.4° (QLC2), 
35.3° (TBM1) and 32.3° (TBM2) in the four scenarios. The deviation of 9 12 from this value, 
A#i2 = 9 12 — 9i2 differs in all the scenarios only through a factor of sin 2#f 2 sin 2 9^ [see 
eq. ( 1481) ]. which is equal to unity within 5% for all the four scenarios. As a result, the 
allowed regions for QLC2 and TBM2 are nearly identical, and larger than those for QLC1 
and TBM1, the last two regions being also almost identical to each other. 

Note that, for all the scenarios, the larger the value of tan/3, the closer the value of 
a 2 needs to be to tt. This statement is true even if we use the laboratory values of mo and 
«2. Moreover, the region with tuq tan/3 > 4.4 eV (i.e. motan/3 > 3.1 eV) is disallowed for 
all values of ct2 and for all scenarios. 
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FIG. 3: Constraints from the 3<r allowed range of 612 on the mo tan (5 (eV) - (X2 (radians) parameter 
space for the four scenarios. The regions above the contours are excluded by data for that particular 



scenario. 



B. 



'13 



to discriminate among scenarios 



The four scenarios that we consider here give different predictions for the values of #f 3 : 
8.9° (QLC1), 1.6° (QLC2), 0° (TBMl) and 3.1° (TBM2). In the absence of RG running, 
therefore, a discrimination between some of these scenarios should be possible in the near 
future. For example, if QLC1 is realised in nature, the value of #13 would be accessible to 
the current generation of experiments. However, the value of #13 changes with RG evolution 
and can either increase or decrease depending on the values of Majorana phases, as can be 
seen from ( )48l) . It is conceivable that the allowed ranges of #13 values for all the scenarios 
will then overlap and the power of discrimination will be lost. It is thus worthwhile to check 
whether one retains this discrimination capability in spite of the RG running. What helps 
in this is the fact that at higher values of Wq tan /3, where one expects large RG effects, not 
all «2 3 values are allowed: the observed values of 612 (as shown in Sec. IV Al) as well as of 623 
restrict the values of the Majorana phases, which in turn restrict the allowed values of #13. 

In Fig. IH we show the pemitted values of #13 in the four scenarios, subject to the con- 
straints of the 3a allowed current ranges of the mixing angles. With the current constraints, 
it should be possible to distinguish between QLC1 and the other scenarios in the next round 
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ml v tan(3 (eV) 



FIG. 4: Values of #13 in the QLC scenarios (left panel) and the TBM scenarios (right panel) as a 
function of mo tan /3, allowed with the current constraints on the mixing angles and mass squared 
differences. The shaded region represents the allowed values of #13 at the laboratory scale in each 
scenario. The thick bands represent upper limits on ttjq tan /3 from current measurements. 



of experiments 
be excluded. 



probing #13. For example, if #13 is shown to be less than 6°, QLC1 will 



The scenario TBM2 may be distinguishable from the remaining two cases if ttiq tan (3 < 2 
eV (i.e. mo tan (3 < 1.4 eV). So this needs, in addition to the information on # 13 , stronger 
constraints on the values of mo or tan f3. A more accurate measurement [40| of 9 12 will also 
indirectly limit the extent of RG evolution and may enable one to distinguish TBM2 from 
the other two scenarios, viz. QLC2 and TBM1. The last two scenarios have very similar 
predictions on # 13 . This is on account of the close values of #^ 3 in the two scenarios, and 
the nearly identical evolution A# 13 = # 13 — 9± 3 in them, whatever the values of , a; 3 , ttIq 
and tan/?. From (jUJ), the value of A#i 3 differs in all the scenarios only through a factor of 
sin 29^2 sin 29£ 3 [see eq. flUDL which is equal to unity within 5% for all the four scenarios. 
Therefore, distinguishing between these two scenarios on the basis of a measurement of 9 13 
alone will be difficult unless the values of the other parameters - mfc, tan/3, , a 3 - are also 
known to a good accuracy. 

Though the constraints on wig tan j3 and in sec. IV Al are logarithmically sensitive to 
the choice of the high scale A, the limits on #13 in a given scenario are almost indepenent of 
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A. This is because the maximum allowed value of |A0j3| is dictated mainly by the maximum 
allowed value of A9\2 through \ A9i 3 \/ A9i 2 ~ | ^13 1 / ^12, whereas the maximum allowed value 
of A6*i2 in turn is an experimentally determined quantity, quite independent of A. 

C. Effect of RG evolution on #23 

The predictions of all the scenarios for 623 are almost identical. The values of 9^ are 
very close: 42.1° (QLC1), 43.4° (QLC2), 45° (TBM1) and 42.7° (TBM2). Moreover, the 
deviations A9 23 = 9 23 — 9 23 are almost independent of the scenario, but depend on the 
value of a 3 . This may be shown as follows: at large values of tUq that are needed to have 
significant RG running, at the high scale \p\\ <C 1, so that T A pa l/p% and the value of a 2 
is restricted to be very close to n. As a result, (J4"5]l gives 

A9 2 3 w -jf sin 20^(1 - cos 20^ 2 cos a^) (56) 
Pa 

which is almost independent of the scenario on account of similar values of 9 23 and 9± 2 , 
but may vary by ~ cos2#f 2 pa 30% depending on the value of a 3 . As a consequence, a 
measurement of #23 cannot discriminate between the four scenarios unless it is accurate to 
the level of a degree. However, if a scenario has already been identified, the measured value 
of 9 23 will restrict the allowed values of the Majorana phase a 3 ■ 

VI. DISCUSSION AND SUMMARY 

The current data on neutrino masses and mixings angles are consistent with symmetry- 
based schemes like quark- lepton complementarity (QLC) or tribimaximal mixing (TBM). 
These scenarios predict specific values of the neutrino mixing angles, which need to be 
compared with their forthcoming more accurately measured values in order to confirm or 
exclude a particular postulated symmetry pattern. However, the symmetry relations need 
to be imposed at a high scale, e.g. the seesaw scale ~ 10 12 GeV, where the neutrino 
masses originate. Radiative corrections to the neutrino masses in general do not respect the 
symmetries involved in QLC or TBM. As a result, predictions of the neutrino parameters 
measured at laboratory energies become different from those given by these symmetries at 
the high scale. It is therefore necessary to obtain low scale predictions of these scenarios. 
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We have calculated radiative corrections to the neutrino mixing angles in the context 
of the minimal supersymmetric standard model (MSSM). We have taken the low scale to 
be the supersymmetry breaking scale ~ 10 3 GeV and have neglected threshold effects, if 
any, during this renormalisation group (RG) evolution. We have presented a technique to 
calculate the deviations A0y of the mixing angles from their high scale values. This technique 
yields analytically transparent results where the errors, caused by the approximations made, 
are small and under control. The analytic treatment clarifies the dependence of the RG 
running of neutrino mixing angles on currently unknown parameters like the mass scale 
mo of neutrinos, the value of tan/3 in MSSM and the values of the Majorana phases. We 
have also solved the RG equations numerically to confirm that the results are indeed closely 
approximated by our analytical expressions. 

We have pointed out certain important patterns in the RG evolution of the mixing angles 
that are valid in any scenario. The RG running of Q\i always increases its value from the 
high to the low scale. Therefore, if a scenario predicts the value of at the high scale 
A, the low scale measurement must be dyi > #f 2 in order for the scenario to stay valid. 
Similarly, the value of 823 increases (decreases), while running from a high to a low scale, 
for a normal (inverted) neutrino mass ordering. The value of 8 13 is controlled not only by 
mass ordering, but also by the values of the Majorana phases, depending on which it may 
increase or decrease with RG running. 

We have considered two versions of the QLC principle (QLC1 and QLC2) and two ver- 
sions of the TBM scheme (TBM1 and TBM2), whose predictions at the high scale are 
consistent with the measured neutrino mixing angles at laboratory energies. The PMNS 
mixing matrices at the high scale predicted within these scenarios are V,l KM U u ^ m (QLC1), 
^bmVcKM (QLC2), Uvtbm (TBM1) and VJ L U u>t hm (TBM2) respectively, where V C km is the 
CKM matrix, U u ^ m the bimaximal mixing matrix, U Utt bm the tribimaximal 12j mixing ma- 
trix, and Vil a charged lepton mixing matrix inspired by the Georgi-Jarlskog relation 36] 
at the GUT scale. 

We summarise our findings in three items: 

(i) The RG running of the mixing angles should not be too large lest the low energy 
values differ too much from their high energy predictions. Since # 12 is the most accurately 
measured angle currently, that puts strong constraints on the allowed values of tan (5 as 
well as on the Majorana phase a^- It is observed that, for m^ tan/?>2 eV (i.e. mo tan /?> 1.4 
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eV) , the allowed range of a 2 is severely restricted. In all the scenarios, Wq tan j3 has to be 
< 4.4 eV (i.e. mo tan (3 has to be < 3.1 eV) for consistency with data. Moreover, the larger 
the value of m tan/3, the closer to n has to be the value of a 2 , and hence that of a 2 at 
the laboratory scale, cf. eq. flo^l) . Thus we find a preference for the approximate equality 
rju ~ — m 2 (especially for large tan/5), as also suggested by considerations of leptogenesis 



2l| . A reduction of errors on the 9\ 2 measurement would decrease the allowed area in the 



ttIq tan p - a 2 plane. Moreover, since 9i 2 always increases from higher to lower scales, if its 
value is measured to be smaller than what is predicted at the high scale in a scenario, that 
particular scenario would get excluded. All the scenarios considered in this paper would be 
excluded if 9\ 2 were measured to be less than 32° . 

(ii) The measurement of #13 is most likely to serve as a discriminator among the four 
scenarios considered here. The predicted values of #13 in these scenarios at the high scale 
are 8.9° (QLC1), 1.6° (QLC2), 0° (TBM1) and 3.1° (TBM2). RG running can modify the 
value of #13 in either direction; however, the restrictions on m^ tan/3 and a 2 from the 9i 2 
measurements limit the extent of this modification. We find, for example, that the value 
of #13 in QLC1 cannot be less than 6° (3cr), whereas in none of the other cases can #13 be 
as large as 6° within 3a. Neutrino experiments during the next decade should be able to 
measure the value of #13 if it is greater than m 5° or to put an upper bound of ~ 5° on it. 
In either case, the scanario QLC1 will be distinguishable from the others. Both QLC2 and 
TBM1 predict almost identical #13 ranges: #13 < 3° (3a). The allowed 3a range of 6*13 for 
TBM2 overlaps with the QLC2/TBM1 range for m^tan/? > 2 eV (i.e. m tan/3 > 1.4 eV). 
Limiting mo tan (3 to lesser values would also help in discriminating between TBM2 on one 
hand, and QLC1/TBM2 on the other. 

(iii) It is not possible for a #23 measurement to discriminate among the four scenarios 
unless it is accurate to the level of a degree. However the value of #23 within any scenario 
is strongly dependent on the Majorana phase ctg . Therefore, if a scenario has already been 
identified, the measured value of #23 will restrict the allowed values of . 

In conclusion, we have shown how the high scale predictions on neutrino mixing an- 
gles get modified with RG running in MSSM for four symmetry-inspired scenarios that are 
consistent with the current neutrino data. With a combination of analytical insights and 
numerical calculations, we show that this limits the allowed ranges of parameters like mg, 
tan p and the Majorana phases. We also indicate the extent to which future measurements 
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can discriminate among various scenarios and how the values of the parameters may be 
further restricted. 
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